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We analyze the static QCD potential in the distance region 0.1 fm < r < 1 fm. We combine 
most recent lattice computations and perturbative computations of the potential, in the frame- 
work of operator-product expansion (OPE). We determine simultaneously the non-perturbative 
contribution to the potential, 8E\js(r), and the relation between the lattice scale (Sommer scale) 
and Ajyjg in the quenched approximation. We find that (1) large part of the short-distance linear 
potential belongs to the perturbative Wilson coefficient, (2) 5£us( r ) — is disfavored, and (3) 
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1. Introduction 

For some years, the static QCD potential VqcdM has been successfully computed in lattice 
simulations. For instance, the distance r of the potential computed in some of most recent lattice 
computations in the quenched approximation ranges from about 0. 1 fm to beyond 2 fm. Computed 
results by different groups, when superimposed with one another, come more or less on a single 
curve with fairly good accuracy, showing very good stability of the lattice predictions (see Fig. [I](a) 
below). In particular, they show that numerically the QCD potential tends to a linear potential 
smoothly at large distances, in accord with the confinement picture of the quarks. 

In this paper, we analyze the potential in the distance region which is smaller than (but not too 
small as compared to) the typical hadron scale Aq^ d . More specifically, we consider the region 
0.1 fm < r < 0.5-1 fm if expressed in physical units. This is the region that is relevant to spec- 
troscopy of the heavy quarkonia, such as bottomonium and charmonium. In this region, it is known 
empirically that the static potential can be approximated well by a Coulomb + linear form. 

In this distance range, accuracy of the perturbative predictions for the static potential improved 
drastically around year 1998 It was proposed to subtract IR renormalon contained in the r- 
independent (constant) part of the potential, e.g. by computing the total energy 2ra po i e + Vqcd(^) 
after reexpressing wi po i e in terms of the MS mass or by computing the force. Then, we find dramatic 
improvement in convergence of the perturbative series, as well as much more stable perturbative 
predictions against scale variation. Moreover, once we obtain much more accurate predictions, 
good agreement with phenomenological potentials and with lattice computations of the potential 
has been observed in the distance range 0. 1 fm < r < 0.5 fm; by now, this has been confirmed by 
several groups [f|,[3],^|]. 

In this article we review our recent work in which we analyze the QCD potential in the 
framework of operator-product expansion (OPE) [Hj. We assemble most recent developments of 
the perturbative predictions and lattice computations. Using OPE, one can separate perturbative 
contributions from non-perturbative contributions in an unambiguous manner. In this way, we are 
able to determine (1) the size of non-perturbative contributions to the potential, and (2) the relation 
between and lattice scale (Sommer scale). In practical applications, the latter determination 
would be particularly important, since it can be used to determine cts(Mz) with high accuracy in 
near future, when lattice simulations incorporating dynamical quarks with realistic masses will 
become available. 

2. OPE of the QCD potential 

The OPE of the static QCD potential was developed around year 1999 within the frame- 
work of potential-NRQCD effective field theory [f7|]. Conceptually it is close to OPEs of other 
physical quantities, with which one may be more familiar. When there is a small parameter as 
compared to the typical hadron size (in our case the distance r <C ^qqd between static quark and 
antiquark), one constructs a series of operators as an expansion in r. Each operator has a Wilson co- 
efficient, which is perturbatively computable, whereas a non-perturbative contribution is contained 
in the matrix element of the operator. The OPE of the QCD potential reads 

V QC D(r)=V s (r) + SE us (r) (2.1) 
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Here, Vs(r) denotes the Wilson coefficient corresponding to the leading UV contribution to 
the potential (whose operator is the identity operator). Vs(r) is referred to as "singlet poten- 
tial," is perturbatively computable, and has a Coulombic behavior (with logarithmic corrections) 
~ — l/|rlogr| at short-distances. We note that Vs(r) is different from the perturbative expansion 
of VQCD( r ); ai l IR contributions (IR divergences and IR renormalons) have been subtracted and 
absorbed into a non-perturbative contribution, so that the perturbative expansion of V$(r) is more 
convergent than that of Vqcd(^). See [Q] for the precise definition of Vs(r). 

On the other hand, SEu$ (r) denotes the leading non-perturbative contribution. It is given as 
an integral over time of an operator matrix element. The operator is local in space but non-local 
in time. The behavior of SE\js(r) at small-distances is known exactly to be £?(AQ CD r 3 ) when r is 
very small (when Vo—Vs& CAOCs/r » Aq CD holds). It is, however, doubious whether this strong 
condition is met in the distance region of our interest, 0.1 fm < r < 1 fm. In this distance region, 
where it is more likely that Vo—Vs~ Aqcd, the behavior of <5£ijs( r ) cannot be predicted model- 
independently. According to some models it is estimated to be 0(Aq CD t 2 ). In any case, in the 
following analysis, we will use only the fact that 5£"us (r) — > as r — > 0, which is true because at 
very small distances 8E\js(r) = ^(AQ CD r 3 ). 

Intuitively one can interpret in the following way. It is the contributions of gluons, whose 
wavelengths are smaller than r, that generates the Coulomb singularity of the potential as r — > 0. 
On the other hand, non-perturbative contributions can be regarded as contributions of gluons whose 
wavelengths are larger than r. Such gluons see only the total charge of the system, hence, as r — > 0, 
they decouple from the color-singlet system. Therefore, their contributions vanish as r — > 0. 

3. Comparison of Vs(r) and recent lattice computations 

In order to make a numerical cross chek of the prediction of OPE, we first compare the per- 
turbative predictions of the leading Wilson coefficient Vs(r) and recent lattice results. In Fig. [T](a) 
we plot lattice data of the QCD potential from three different groups @ g|. In the same figure, the 
perturbative predictions of V$(r) are plotted up to the next-to-next-to-next- to-leading logarithmic 
order (NNNLL). As a salient feature, we see that the perturbative predictions follow the lattice data 
up to larger distances as we increase the order. 

Following input parameters were chosen in depicting this figure. First, we note that in order 
to make a comparison between Vs(r) and lattice data, the relation between A^j§- and lattice scale 
(we take the Sommer scale ro) is needed. This is because, a priori each lattice data set is given in 
units of ro, whereas perturbative predictions are given in units of A^. Thus, conversion of units 
is needed to compare them in common units. Here, only in this section, we use the central value 
of the conventionally known relation r A^ p = 0.602 ± 0.048 Other inputs are: a s (Q) = 0.2, 
ni = 0; at NNNLL, the 3-loop non-logarithmic term is not known yet, only some estimates exist, 
hence, we used Pineda's estimate for [Q] and varied it within a range including other estimates 
as a part of our estimates of uncertainties; we chose C+L scheme in computing Vj(r). 
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To verify the prediction of OPE, a closer inspection is needed. Data points in Fig. [j](b) show 
the differences of Vs(r) and the lattice data [||], where the vertical scale is magnified as compared 
to Fig. |l](a). At LL, one can still see a Coulombic contribution in the difference, but as the order 
is increased, a Coulombic contribution vanishes and they tend to be regular at the origin. Solid 
lines in the same figure show fits of the data points by cubic polynomials, using data points at 
rA^ p < 0.5. One sees that the solid lines approximate the data points up to larger distances as 
the order increases, which indicate that these data points become regular around the origin as the 
order increases. These features support the OPE prediction, that the non-perturbative contribution 
8E\js(r) vanishes as r — > 0. We note that such behavior can be observed only when we take the 
difference between Vs(r) and the high quality lattice data, since in Fig. |j](a), perturbative predictions 
become so steep toward the origin that any steep curve tends to go through the lattice data points. 



4. Determination of 8Eu S (r) and r A^ 

In the light of the observation in the previous section, we perform a simultaneous fit to de- 
termine the non-perturbative contribution 5£\js ( r ) an d the relation between lattice scale and A^-, 
x = roA^ p . (In the previous section we chose a specific value for x determined from another 
source, but in this section we determine its value using the QCD potential alone.) 

Let us explain why we have a high sensitivity to the value of x. It constitutes the essence of 
our analysis in this article. We compare the lattice data for the potential and perturbative prediciton 
of Vs(r) in some common units. Here, we choose to compare in units of A^, hence we convert the 
lattice data into these units by using any chosen value of x = r$ A^°| p . Let us write 

VMt(r;x)-V s (r) = [V Mt (r,x tme ) -V s (r)} + [V Mt (r,x) -V latt (r,x tme )]. (4.1) 
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Figure 2: Allowed regions in the parameter space (Ai,A2) at 68% and 95% CL, when Viatt^;.*) — Vs(r) is 



It is a trivial equality, where x tme denotes the true value of x. The first bracket [...] on the right-hand- 
side coincides with the non-perturbative contribution 8E\j$(r), hence it vanishes as r — > 0. On the 
other hand, the second bracket on the right-hand-side takes a "Coulomb"+linear form as long as 
x 7^x trae . (Here, "Coulomb" includes logarithmic corrections at short-distances.) It follows from 
the fact that Vi at t(r;x) takes a "Coulomb"+linear form and that [Viatt^x) — VJatt^JCtrue)] i s merely 
the difference of these functions after rescaling Vi att (r,x) in the vertical and horizontal directions. 
Hence, the second term contains a Coulombic singularity as r — > unless x = x tme . It means that 
in order to find the true value of x, we should adjust x such that the leading Coulombic behavior 
vanishes in V\ an (r;x) — Vs(r) as r — > 0. That the fit is determined by the leading singular behavior 
guarantees a high sensitivity to x. 

As for 8Eus(r), we simply fit it by a quadratic polynomial using the data points at 7"A^ P < 0.5, 
since 8E\j$(r) is regular as r — > 0. 1 

We show the results of the fit for 8E\js(r) in Fig. ||[ The central value of the coefficient 
of the linear potential is about 1 in units of (A^ P ) 2 - We note that the coefficient of the linear 
potential (string tension) determined from the large distance behavior of lattice results is about 3.8 
in the same units. Hence, our result (central value) indicate that the non-perturbative contribution 
contains only about one quarter of the string tension, as a component of the linear potential in the 
short-distance region. Three quarters of the string tension belongs to the perturbative contribution 
Vs(r). The present error is somewhat large, so it is still consistent at 95% confidence level that 
either all of the string tension is contained in Vs(r) or only 1/3 of it is contained in Vs(r). 

Another notable feature is that (in the schemes which we examined) 8E\js(r) = is strongly 
disfavored. It is seen from Fig. ||, where the origin is excluded from the allowed region. It means 
that there is a non-vanishing non-perturbative contribution. In fact, this is a new result, since in 
a previous similar study by Pineda [0], 8E\js(r) contained a larger error and was still consistent 
with zero. The major difference of the analyses, which led to the improvement of the bound in our 
analysis, is that we determined x simultaneously from the fit. 

We also obtained 

* = 0.574 ±0.042. (4.2) 
'We also tried a cubic fit, where the results of the fit did not change considerably. 
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It should be compared with the conventional values determined by the Schrodinger functional 
method: 

x = 0.602 ±0.048 [§, (4.3) 
x = 0.586 ±0.048 [§. (4.4) 

The errors are of the same order of magnitude, and the values are consistent with one another with 
respect to the errors. We emphasize that our determination uses a completely independent method, 
and that the mutual consistency between these values is quite non-trivial. At the moment, the error 



of our x in eq. (4.2) is dominated by errors of the lattice data of the potential. 



5. Conclusions 

We analyzed the static QCD potential in the distance region 0. 1 fm < r < 1 fm. We combined 
the perturbative predictions up to NNNLL and the most accurate lattice computations, in the frame- 
work of OPE. In this way, we determined the non-perturbative contribution 8E\js(r) and roA^ p , 
for ni =0 (quenched approximation). Our conclusions are as follows: 

• Most of the linear potential at r < 1 fm (~ 3.8 A^— r) is included in perturbative prediction 
of the Wilson coefficient Vs(r). This is a scheme-independent result. 2 

• 8E\js(r) = is disfavored. This is a scheme-dependent result, but it holds in different 
schemes which we examined (which we consider to be reasonable schemes). 

• roA^ p = 0.574 ± 0.042. It is consistent with the conventional values and the error is of 

u MS 

comparable size. This provides a new method for its determination. 
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2 Scheme dependence enters at Gij 1 ), hence it does not affect the linear potential. 
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